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1. Introduction 

General relativity has modified our understanding of the physical world in depth and 
has altered some among the most fundamental notions we use to describe it. During the 
last ten years, the effort to understand the combined consequences of this conceptual 
revolution and quantum mechanics has lead to loop quantum gravity, a predictive 
quantum theory of the gravitational field, whose theoretical results can be, in principle, 
empirically tested 0. There are other areas in our understanding of nature, however, 
where the consequences of the general relativistic conceptual revolution have not been 
fully explored yet. Among these is statistical mechanics. To be precise, thermodynamics 
and statistical mechanics on a fixed curved spacetime have been much studied (see, for 
instance, 0); but not much is known on the possibility of developing thermodynamics 
and statistical mechanics of a fully general covariant system, in particular, a system 
including the gravitational field. 

Here, we begin to address this issue. Specifically, we study the following problem. 
Consider a simple physical system, s, with a finite number D of degrees of freedom. 
Assume that s is described by a fully constrained Hamiltonian system. That is, its 
dynamics is not given by a Hamiltonian, but rather by M first class constraints. 
Physically, this means that we do not understand the dynamics of s in terms of 
the evolution of D dependent Lagrangian variables (or 2D phase space variables) 
as functions of a single preferred independent external time variable t; rather, we 
understand the dynamics as the relative evolution of 2(D+M) phase space variables with 
respect to one another - all the variables being on the same footing. The dynamics fixes 
relations between these variables, so that by knowing some of them we can predict the 
others. Such a simple system encodes an critical feature of general relativistic systems: 
the absence of a preferred time variable, and the relational aspect of evolution. Now, 
consider a macroscopic system S composed by a large number of component systems, 
each one identical to s, and interacting weakly. Can we use statistical mechanics to 
describe macroscopic properties of S? Notice that there is no time variable in the 
description of S, therefore no notion of thermalization 'in time'; there is also no notion of 
energy, and thus no obvious way to define a canonical or microcanonical ensemble. If we 
arbitrarily choose one variable in S as the physical time (that is, if we 'deparameterize' 
the system), and then use conventional statistical techniques, our results are going to 
depend on the choice of time, and therefore to be possibly unphysical. Is there anything 
we can nevertheless say, about the macroscopic behavior of this system? Can we still 
apply thermodynamical or statistical mechanical techniques? 

These questions are relevant in a strong- field gravitational context, whenever a 
preferred time and a conserved energy are not defined. Of course, if we consider a system 
with a notion of time and with conserved energy, we expect that temperature and energy 
will recover their traditional role. This is the case, for instance, of an asymptotically 
flat gravitational field; in this case the Hamiltonian is given by suitable boundary terms 
and the observables at infinity evolve in the Lorentz time of the asymptotic metric. The 
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general theory we present here will have to yield standard results in this case. However, 
what about the situations in which there is no conserved energy and no preferred time? 
For instance, as far as we know, our universe might very well not be asymptotically flat. 
Alternatively, we may be interested in a system with a strong (dynamical) gravitational 
field, and have no access to an external asymptotic region. In particular, consider a "high 
temperature" early-universe regime. This is usually described in terms of fluctuations 
around a background metric; is there a genuinely general covariant description of this 
physics? And what is temperature in this context, if we do not fix a background metric? 
Certainly, it is difficult to even define what statistical mechanics is if we do not have some 
notion of energy conservation; but does this mean that in all gravitational systems in 
which there is no conserved energy (most of them!), we have to renounce using statistical 
methods? 

These questions have not yet been addressed in the literature, as far as we know. 
An attempt to study certain aspects of the foundations of general covariant statistical 
theory is in Refs. [|| |3j]. In these works, the question addressed is whether a preferred 
time flow, having the thermodynamical properties that we ascribe to physical time, can 
be derived from the statistical mechanics of a covariant system. The answer is positive, 
and the flow turns out to be dependent on the statistical state. The relation flow/state 
reflects a very general operator algebra structure (Tomita-Takesaki theorem), and raises 
intriguing physical issues, in particular in view of powerful mathematical uniqueness 
results about the flow (Connes' Cocycle Radon-Nikodym theorem). Here, on the other 
hand, we are not concerned with the emergence of a time flow. Instead, we address 
directly the issue of a statistical description independent from any notion of time. 

Furthermore, Refs. |3|, |J take Gibbs' || (and Einstein's ||!) point of view on 
statistical mechanics: a statistical state is described by a distribution over the phase 
space T of the composite system S (in Ehrenfest's terminology, over the T-space 0). 
The state represents the distribution of S"s microstates over many imaginary copies of 
the system, all in the same macrostate. Here, on the contrary, we use Boltzmann's 
original point of view ||: we assume that S is composed by a large number of identical 
subsystems s. The statistical state is then described by a distribution over the phase 
space 7 of the component system s (over the ^x-space, in Ehrenfest's terminology). This 
gives, for each state of s, the expected number of component systems that are in that 
state. 

Of course, we do not expect any of the well known subtleties and conceptual 
difficulties of statistical mechanics to be solved by applying it to covariant systems. Here 
we are not concerned with the old problems in the foundations of statistical mechanics, 
but only with the specific new problems -and new beauties- that emerge in trying to 
extend the general relativistic revolution to statistical physics. 

Our main result is the following. We argue that, under appropriate conditions, 
the statistical mechanics of a system S composed by many constrained systems s 
is well defined. In particular, statistical mechanics is not necessarily tied to the 
concept of energy, or to a preferred time flow. Accordingly, general covariant statistical 



4 



mechanics is not governed by the notion of temperature. Instead, intensive macroscopic 
parameters are determined by the properties of Boltzmann's thermalizing interaction 
among the individual component systems. In the course of the paper, we develop the 
basis of covariant quantum statistical mechanics and define the intensive and extensive 
thermodynamical quantities. 

We begin by recalling the properties and the physical interpretation of the 
parameterized systems in Section 2. We then give the main discussion on the foundations 
of covariant statistical mechanics in Section 3, and a simple example in Section 4. We 
discuss the statistical mechanics of a gas of free relativistic particles in Section 5, and 
we comment and summarize in Section 6. 



2. Presymplectic systems 

We consider fully constrained systems, with a finite number of degrees of freedom, and 
with first class constraints | i~0| . Their dynamics is obtained from the action 



S[q\ Pi , X m ] = Jdr \j^Pi- A m C m (<f,^)} , (1) 

which is invariant under arbitrary reparametrizations of the parameter r. The parameter 
r is unphysical and unobservable, like the time coordinate in general relativity. The 
unreduced, or extended phase space j ex is coordinatized by the canonical pairs (q t ,Pi); 
i = 1, 2, N. The canonical 2-form on j ex is uj ex = dpit\dq % . The pair (T ex , u ex ) forms a 
symplectic space. The variation of the action with respect to the canonical coordinates 
q l ,Pi gives the equations of motion 

H_ = } m dC m{.q\Pi) 

dr dpi 

dpi dC m {q\pi) 

while the variation of the action with respect to the Lagrange multipliers A m gives the 
constraint equations 

C m =C m tf, Pi ) = 0, ro = l,2,..,M. (3) 

Thus, the dynamics of the system with respect to r is the unfolding of the gauge 
symmetry generated by the first class constraints, i.e., dynamics is gauge. 
The first class constraints satisfy, in general, a "non-Lie" algebra 

{C m ,C n } = C mn \q i ,p l ) Q, (4) 

and the number of independent physical degrees of freedom of the theory is D = N — M. 
The constraint surface 7 in 7 ex defined by the constraint equations @ is a (2D + M)- 
dimensional manifold. The restriction u of u ex to the constraint surface 7 is of rank 
2D. The M null directions of u are the infinitesimal transformations generated by the 
constraints. They define the gauge orbits on 7. The physical phase space 7^ is the space 
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of these orbits. This is the space of the physically distinct solutions of the equations of 
motion. 

The space (7, u) is a presymplectic space, which contains the full dynamical 
information about the system. Hence dynamical systems in this form are also called 
'presymplectic systems'. 7 can be parameterized by the set of independent coordinates 
{q a ,p a ,t m ), where (q a ,p a ), a = 1,2, ...,D are canonical variables that coordinatize the 
physical phase space 7 p h, and t m , m = 1,2,..., M coordinatize the orbits. In general 
this coordinatization can hold only locally, and different charts may be needed to cover 
the entire space. 

Any conventional dynamical system with phase space (7^, u p h = dp a A dq a ), and 
Hamiltonian H = H(p a , q a ) can be represented as a presymplectic system as 



where t is the coordinate in R, and corresponds to the external time variable. The 
difference between the conventional formulation and the presymplectic formulation is 
only in the fact that this time variable is treated on the same footing as the other 
variables. As a concrete example, we may imagine that H is the harmonic oscillator 
Hamiltonian describing the small oscillations of a pendulum, while t is the reading of 
a physical clock. Then the presymplectic system (|5]) describes how two equal-footing 
physical variables (the pendulum amplitude and the clock reading) evolve with respect 
to one another. In general covariant systems, such as any general relativistic system, 
this 'equal footing' status between all physical variables is an essential feature of the 
theory. It expresses the major physical discovery of general relativity: the complete 
relativity of spacetime localization. 

Note that the canonical coordinates q a , and p a are the physical observables of the 
system. They are gauge-invariant. They satisfy {q a ,pb} = 8% on the physical phase 
space. In these coordinates, the physical symplectic form on 7^ is uo p h = dp a A dq a . The 
general solution of the equations of motion is simply given by the embedding equations 
of the orbits in j ex , that is 



Each set {q a ,p a ) determines a solution; along each solution, the quantities (q l ,p%) depend 
on the M parameters t m (instead than just on a single time variable) because of the 
gauge freedom in the evolution. The inverse relations of (|6])-(|7|) give the dependence of 
the physical observables q a , and p a from the original coordinates 



(7 = l P h* R, w = ^ P h ~ H(p a , q a ) A dt) 



(5) 



q l = q\t m ;q a ,p a ), 
Pi=Pi(t m ;q a ,p a ) . 



(6) 
(7) 



q a = q a (q\Pi) 

Pa =Pa(q\Pi) 



(8) 
(9) 



as well as the orbit coordinates t 



ni 



t m = i m {q\p l ) 



(10) 
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The quantities (§) and commute with all the constraints, and provide a complete set 
(in the sense of Dirac) of gauge-invariant observables. Every other physical observable 
can be obtained from them. 

Let us recall how evolution can be obtained from the basic observables (H) and (|9|) 
TT| , |12f . If we plug the gauge variables ( TTOj ) into the full solution (|(J) and (|7|) we obtain 



the equations 

q l = q\t m (q\p t );q a ,Pa), (11) 
p l =p l (t m (q\p l );q a ,Pa)- (12) 

In general, 2N — M of these equations are independent. For each physical state of the 
system, determined by the value of (q a ,p a ), these equations define an M dimensional 
subspace in the phase space. Therefore each state determines a set of relations on the 
original phase space variables. These relations represent the dynamical information on 
the system; they provide the full solution of the dynamics in a gauge-invariant fashion 

0- 

In particular, we might arbitrarily choose a set of M coordinates q m (or momenta 
p m ; or a combination of both) as independent 'clock and position' variables, and express 
the evolution of the remaining set of coordinates and momenta as functions of these 
q m for any physical state (q a ,p a )- For each fixed numerical value q m of the coordinates 
q m , we have a well defined gauge-invariant observables in 7^. For instance, let us 
chose (arbitrarily) q 1 as a dependent 'partial' observable[|, and the next M of the q l, s, 
as independent 'partial' observables, or 'clock and position variables'. That is, let us 
choose m = 2,...,M+l. Pick M fixed numerical values q m for the M variables q m . 
Generically, this fixes uniquely a point on every orbit. The value Qi m of q 1 on this point 



depends on the orbit, and can be obtained from (|8|-|iOD. Let it be 

Q\m = Q\ m {q\p a ). (13) 

Here all the q l are partial observables, while Q^ m is a complete observable. The function 
([jjD is gauge- invariant, well defined on 7^ and expresses the relative evolution of q 1 , as 
a function of the q m , m = 2, . . . , M. It is called an 'evolving constant of the motion', or 



simply a 'relational observable' [11 



The quantum theory can be constructed by imposing the quantum constraints 
on the unconstrained Hilbert space 7i (or some suitable extension of the same if 
the constraints have continuum spectrum). The space of solutions of the constraint 
equations is the physical Hilbert space TC p h ys of the theory. (If H p h ys is not a subspace 
of H a scalar product is determined in TC p h ys by the requirement that the self-adjoint 
observables in 7i which are well defined on 7i p hy S be still self-adjoint.) Generically, 
we expect that out of the operators corresponding to the set of 2D gauge invariant 

§ A 'partial' observable is a physical quantity to which we associate a number, such as time t, position 
x or electric field E. A 'complete observable' is a physical quantity that can be predicted if the state 
is known, or, equivalently, that gives us information on the state, for instance the value E(t, x) of the 
electric field in a certain point a; at a certain time t. For the notions of partial observable and complete 
observable, see fl3|. 
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observables (q a ,p a ), we can define D = N — M commuting operators, a , a = 1,2, D 
forming a complete Dirac set. Assuming for simplicity these have discrete spectrum, 
a basis of physical states is labeled by their quantum numbers n a , a = 1,2, ...,D. A 
general physical state is killed by all the constraints C m \tf)) = |0), and can be written 

as 

W) = ^2 c m,-,n D \nu...,n D ). (14) 

m,...,ri£> 

Physical evolution is described by (Heisenberg) operators corresponding to relational 
classical quantities such as (|I3D . In constructing these operators, ordering and 
consistency problem might, in general, be serious. 

3. Covariant statistical mechanics 

Can we use statistical mechanics methods in a covariant, presymplectic framework? 
Energy plays an important role in statistical mechanics, and here there is no energy. 
Statistical mechanics relies on the idea that systems thermalize to equilibrium in time. 
What is thermalization in a covariant context, in which there is no external time 
variable? To address these questions, our strategy will be to recall Boltzmann's logic, to 
rephrase it in the language of the presymplectic formulation of a conventional system, 
and from here, to extend it to presymplectic systems that do not correspond to a 
conventional system. 

Consider a Boltzmann gas in a closed box. The gas is composed by a very large 
number M of identical molecules. Begin by considering each molecule as free. Let 
7 be the phase space of a single molecule. For instance, if we neglect rotational and 
vibrational motion, we may assume 7 to be six dimensional. Since the molecule is 
assumed to be free, its motion is very simply described by a free motion in 7. The phase 
space T of the entire gas has dimension QKf . The motion of the entire gas is described 
by a simple motion in V as well. Under these assumptions, the gas does not thermalize, 
and we cannot use statistical methods. For instance if we started with all the molecules 
bouncing up and down within the right half of the box, they would continue to do so 
forever, never expanding to the left-hand part of the box. To have thermal behavior, we 
need the particles to interact. However, taking the actual physical interaction among 
the molecules into account complicates the dynamical problem dramatically, and puts 
it far outside our theoretical capabilities. 

Boltzmann's genius found a way in between, by postulating a 'small,' 'thermalizing' 
interaction among the molecules. The molecules bounce, attract and repel in a non- 
trivial manner. In the theoretical description, we simply assume that each molecule is 
still free most of the time, but, once in a while, it interacts with another molecule. We 
are not concerned with the details of this interaction, except for the assumption that 
the interaction is maximally thermalizing, that is, it conserves a minimal number of 
physical quantities. Under this assumption, motion in Y becomes ergodic and we have 
thermalization. As time goes on, the state of the gas will fill up all allowed regions in V. 
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Of course, there are quantities that must be conserved in any interactions, due to the 
homogeneity properties of the spacetime in which the gas lives, such as momentum and 
energy. The presence of the box walls forces the total momentum to be zero, and the 
only non-trivial conserved quantity is the total energy. Anything else is washed away by 
the thermalizing interaction. We assume that time averages are the same as ensemble 
averages, and that under the action of the thermalizing interaction all microstates of 
the gas become equiprobable, with the only constraint given by the value of the total 
energy. Thus macroscopic (microcanonical) states can be labeled by a single parameter, 
their total energy. As is well known, the quantitative consequences of this very delicate 
argument, considered borderline fantasy by Boltzmann's contemporaries, are strikingly 
accurate in a truly impressive range of physical contexts. 

In the course of the dynamics, the motion of a single molecule can be followed 
within its phase space 7. This motion is free for most of the time, but at certain times 
it gets suddenly altered: when the molecule interacts with another molecule. Assuming 
equiprobability, a simple calculation shows then that the time averaged distribution 
of the states of a single molecule, and thus the distribution of the molecules over the 
states, is given by p ~ e~^ H , where H is the free Hamiltonian of the particle and f3, the 
(inverse) temperature, can be computed from the total energy. 

Let us now describe the same system in the presymplectic framework. First, let the 
particles be free. The key difference with the previous description is that a point in the 
physical phase space 7^ does not represent anymore the state of the particle at some 
time. Rather, it represents a single full solution of the equations of motion. (It is like 
a classical analog of a Heisenberg state, versus a Schrodinger state.) Thus, the particle 
motion is now described by a single, non moving, point in j p h, which represents a full 
orbit in 7. Similarly, the motion of the entire gas is given by a single non-moving point 
in T ph , corresponding to a full gauge orbit in T. As there is no time, there is no time 
for moving around. 

However, the magic, once again, happens when we turn Boltzmann's 'small' 
interaction on. The dynamics of the full system is still given by a single non-moving 
point in F p h, or, equivalently, by a single orbit in T. However, what about the dynamics 
of a single molecule? Since the phase space j p h is defined by the dynamics of the system, 
and not just its kinematic as in the conventional case, it seems that the motion of a 
single molecule cannot be described in j ph at all in the interacting situation, because 
j p h is the space of the free motions of the molecule. It seems to be a core difficulty 
However, there must be a way out, since, after all, we are describing the same physics as 
before. Indeed, the way out is provided precisely by the assumptions about Boltzmann's 
thermalizing interaction. Observe that the orbits in T do correspond to free motions of 
the single particles, interrupted by interactions. Each such orbit gives, for every particle, 
a collection of free motions, namely a collection of points in j p h. In other words, what 
the interaction does is simply make the (timeless) state in j ph diffused. A full orbit in 
T p h determines a distribution of points in 7^. Under our assumptions, the density of 
these points must clearly be given by p ~ e~^ H \ 
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What conclusion can we draw from this exercise of re-expressing Boltzmann's ideas 
in a timeless language? The first conclusion is that we can still think in terms of the 
Boltzmann's distribution on the states of the subsystem, even in a timeless context. It 
is true that nothing moves in the physical phase space of a fully parameterized system, 
and so it seems that nothing can ever thermalize. But the effect of the interaction 
between the subsystems can be represented precisely by a distribution on the space of 
timeless non-interacting states§|. 

The second conclusion regards the energy. Why does the energy still play a role, 
when the presymplectic formalism treats the time variable, and thus the energy, just 
as one among other variables? The answer, from the above discussion, is not that the 
energy has any special importance by itself. Rather, it is that we have simply fed into 
the formalism the information that the small interaction between the subsystems washes 
away everything excepts energy. But there is nothing sacred about energy conservation. 
Energy conservation is just a consequence of invariance under time shift, which, in turn, 
is a feature of the homogeneity of the Minkowski solution under time shifts. We have 
learned from general relativity that the Minkowski gravitational field is just one among 
many possible fields. There is no fundamental energy conservation in nature. 

On the other hand, the discussion above leads us to see precisely under which 
conditions we can still use Boltzmann's statistical mechanics in a covariant context. We 
can, anytime we have a system S that can be seen as composed by a large number of 
identical subsystems s, whose dynamics is given by a free part which we understand 
well, plus a 'small interaction' that can thermalize the macrosystem, and conserves, say, 
only some global quantities 0\. We can formalize our conclusions as follows. 

Let (7, oj) be the presymplectic space describing s. Let S, with presymplectic space 
(r, Q) be composed by a large number Af of systems S( n ), n — 1, . . . , Af, all identical to 
s, having presymplectic spaces (7(n), C0( n ))- By this we mean 

r = x„ 7 (n) , (15) 

= y^(n) + U int , (16) 
n 

where x indicates the Cartesian product and Ui nt gives the interaction between the 
subsystems. Next, let us assume that there are L quantities Oi, defined on 7^ (and 
thus on 7) such that the corresponding global quantities 

O l = J2ol n) , l = l,...,L<D, (17) 

n 

t) In some sense, we are dealing here with a 'covariant ideal gas'. As in standard ideal gases, there is 
no specific interaction in the formulae involved. Nevertheless, the ideal gas thermalizes thanks to the 
Boltzmann thermalizing interaction, which is implicitly assumed. If specific interaction terms among 
the single components s were allowed, then it would be interesting to study what thermalization could 
mean in that context, as well as its relationship with the many-body forces required to get separability 
of the whole system S (the cluster decomposition property) Q. This issue is not addressed here and 
we leave it for future developments. 
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are invariant along the orbits of f2, that is 

X(O l ) = 0, l = l,...,L<D, (18) 
for all vector fields X in T such that 

n(x) = o. (19) 

We call these quantities 'conserved'. Finally, we assume that Ui nt suitably thermalizes 
all other variables besides the CVs. This means, precisely as above, that all allowable 
(combined) states of the S( n ) systems are, on average, equally covered, in moving along 
a generic orbit in F. 

Under these conditions, we can straightforwardly construct a covariant statistical 
formalism. A state in T p h, determines, for each sr n \ a distribution p on 7 P h, which 
gives the distribution of 'initial states' of the component system as we move around the 
corresponding orbit of the interacting composite system. For a generic state in T p h, this 
distribution can be computed using conventional statistical techniques, in particular, 
by assuming that the distribution is the one that maximizes the number of possible 
microstates compatible with the given macrostate. The result is straightforward: the 
(unnormalized) probability distribution on the phase space is 

p = e^°\ (20) 

where j l are the intensive thermodynamical parameters that determine the equilibrium 
of the members of the ensemble with respect to the transfer of the quantities [0{\. 

Instead of detailing the classical theory, we discuss directly the quantum theory. 



We use von Neumann's density operator formalism |]14| . We ask the ensemble to satisfy 
a maximum entropy principle. In other words, we ask the quantum statistical entropy 
S per constituent member of the ensemble given by 

S = -k Tr (pmp) , (21) 

to be a maximum under the constraints 

Tr (pO { ) = O h 

Tr p = 1 , (22) 

where 0\ are the quantum operator corresponding to the conserved quantities 0\ and 0\ 
are fixed average values, p is the density operator. The density operator p that fulfills 
these requirements is 

p = Z~ x e"^ 5 ', (23) 

with 

Z = Tr e- y5i , (24) 

the partition function. The thermodynamical parameters 7' can be obtained from the 
conditions ( |2"2"D provided that the matrix have non- vanishing determinant. Clearly, 
they are the parameters that measure the equilibrium of the members of the ensemble 
with respect to the transfer of the quantities [Oj\. 
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In the conventional case of a non-covariant system formulated in covariant terms, 
only one non-trivial quantity is conserved, the energy, and we obtain the standard 
results. In particular, we can consider a gravitational system with an asymptotically 
flat gravitational field; in this case the Hamiltonian is given by suitable boundary 
terms and the observables at infinity evolve in the Lorentz time of the asymptotic 
metric. The system is Lorentz invariant for the asymptotic Lorentz transformations, 
and therefore, in particular, invariant for time translations. Therefore all interactions 
conserve the asymptotic Lorentz energy and the theory considered here reduces to the 
standard results. A specific example of this is given by all the literature on black hole 
thermodynamics an statistical mechanics, in which, in general, the gravitational field is 
assumed to be asymptotically flat. 



4. An example 

As a simple example, we take as component system s a model with two non-commuting 
Hamiltonian constraints and one physical degree of freedom which was studied in ||15|| . 
This model mimics the constraint structure of general relativity. We refer to [15| for all 
details. 

The model we consider is defined by the action 

S[u, v,N,M,X] = ^Jdt[N {Vu 2 + v 2 )+M {Vv 2 + u 2 ) ] , (25) 

where 

Vu=-{U-\u), Vv=^(U+\v); (26) 

the two Lagrangian dynamical variables u = (u 1 ^ 2 ) and v = (v l ,v 2 ) are two- 
dimensional real vectors; N, M and A are Lagrange multipliers. The squares are taken 
in R 2 : u 2 = u ■ u = {u 1 ) 2 + (w 2 ) 2 . The action can be put in the form ([]]), 

S[u,v,p,n,\ m ] = J dr ^u-p+'v^-\ m C m {q\ Pi )] . (27) 

The canonical variables (u, v, p, n) define the eight dimensional extended phase space 
7 ex , with symplectic form u ex = dp A du + dn A dv, also A 1 = N, A 2 = M, and A 3 = A. 
The constraints 



C 1 = 1 -(f-v 2 ), 
1 

■"2 — 

- 2 



C 2 = l(n 2 -u 2 ) 



C 3 = u-p-v-n, (28) 

are first class and define a five dimensional constraint surface 7 in j ex , and uj = u> ex \^. 

A complete set of gauge-invariant quantities is given by the two continuous 
quantities J G R + , 4> G Si and two discrete quantity e, e' = ±1, defined by 

12 12 
up — p u 



| M lp2 _ pl u 2| ' 



12 



12 12 
, TT V — V TT 

| 7r l t ,2 _ ylyj-21 



J = \u p 2 — p 1 !^] , 

12 12 

U V — p TT , . 

= arctan — — — - , (29) 

U V — p L TT L 

These can be taken as coordinates of the physical gauge-invariant phase space. The 
quantity J resembles an angular momentum, and thus it is called as such. 

Let us now consider a large number of systems of this kind which are interacting 
weakly. The composite system dynamics is given by the presymplectic system ([IB]). Let 
us assume, as an example, that u int is a sum of binary interactions in which the sum of 
the two angular momenta, while all other quantities are thermalized. This defines the 
statistical mechanics of the composite system. 

In the quantum theory, we take a complete Dirac set of commuting operators J, e 
and e'. Their spectrum, worked out in is 



J\m, e 


,e') N = 




e,e') N = mh\m,e,e') N , 




e|m,e, 


e% = 


e m, e, 


e') N , 




e'\m, e 


e') N = 


e'\m, e 


,eV 


(30) 


m is a positive integ 


er, and 


e and 


e' take values 1 and -1. 


The states m, e, e') n form a 



normalized basis in the physical Hilbert space of the theory. 



N 



(m, e, e'\rh, e, e) N = 5 mj rh'S e ,eSe',e' ■ (31) 



The quantity J is represented by the operator J, which is a kind of angular momentum 
T5| as we have mentioned. If we assume that the thermalizing interaction conserves the 



total value of J, we have immediately the density operator 

P =L Z~- ( 32 ) 
Using Trp = 1, we get the partition function 

4 

eT R - 1 ' 



m=l 

with u m = 4 because there are 4 states for a given m, and also e~ lh < 1 has been 
used (i.e., positive 'temperature' 7 has been assumed). The angular momentum per 
constituent is given by 

and the entropy per constituent 

S = ka = —kTi (pin p) = k In Z + k^L . (35) 

The parameter 7 characterizes the equilibrium state of the system with the reservoir 
and plays here the role of a temperature. If we had an empirical thermodynamics 
of this system, we could identify this parameter with an empirically determined 
thermodynamical quantity. 
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5. Gas of free relativistic particles 

Before concluding, we discuss a simple case, in which some of the ideas presented 
above can play a role: the case of a gas of free relativistic particles. This is really 
an oversimplified situation, which can be treated with simpler tools; but the illustration 
of this case may be instructive, and can be seen as a check that the theory described 
here is in agreement with other methods in the cases in which other methods can be 
applied. For simplicity, we remain here in the classical context. 

Consider thus a gas of relativistic particles. Can we associate a temperature to 
this gas? What is the statistical state describing these particles? The Hamiltonian 
description of a single relativistic particle can be formulated in a manifestly Lorentz 
covariant fashion as follows. The phase space is coordinatized by the coordinates x M 
and their conjugate momenta p M -that is, the symplectic form is u ext = dx 11 A dp^- and 
the dynamics is given by the constraint C = p 2 — m 2 . The seven dimensional constraint 
surface 7 defined by C = with its induced restriction u of u ext form the presymplectic 
space (7, u>) describing the dynamics of the particle. 

Now, one may say that for this system we know that the time is t = x°, and 
the energy is E = po- Therefore we can apply standard statistical mechanics with no 
difficulties. However, there are two distinct problems. The first is that the entire 
dynamics of the system is contained in the geometry of 7, which has no a priori 
specification of which variable is time and which variable is energy. Thus, can we 
do thermodynamics just on the basis of the actual dynamical laws, without specifying 
which one is the time parameter? However, there is also a second problem, much more 
concrete. Suppose we say that a; is the time variable, and po is the energy. You, on 
the other hand, use a different Lorentz reference frame, and therefore for you the time 
is x'° = A°x^ and the energy is p' Q = AqP m , where A is a Lorentz transformation. If I 
write a Boltzmann statistical state using my definition of energy, and you in yours, do 
we define the same statistical state? It is easy to see that the answer is no, because 



whatever is f3' . So, which one is the correct equilibrium state? 

Let us address both problems in terms of the general theory developed above. The 
key point is that if the gas is formed by particles that are really free, they will never 
thermalize. Some thermalizing interaction is needed in order to reach an equilibrium 
state. Thus, we need some additional physical input (this is the key point). On physical 
grounds, we may for instance observe that our gas of relativistic particles thermalizes 
by means of relativistic elastic scattering. Therefore, the dynamics of a single particle 
is not really free: the presymplectic space describing the dynamics of the system is 
the cartesian product of the spaces of the particles, the total presymplectic form is the 
sum of the individual presymplectic forms plus the interaction term, as in Eqs.([T5] and 
|TB]) . What are the conserved quantities Oi, in the sense of Section III? Namely, what 
are the quantities that are exchanged, but whose total value is conserved, in such an 
interaction? Clearly, they are the momenta p M . Therefore, according to the general 




(36) 
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theory of Section III, the (unnormalized) probability distribution of the states of the 
single particle is 



where the quantities 7 M are the intensive parameters describing the system. 
Straightforward application of standard statistical techniques tells us then that the 
average 4- momentum is proportional to 7^: 



Therefore, 7 M must be timelike, and therefore there is a preferred Lorentz frame in which 
7 = and 70 = (3. This is the frame in which the center of mass of the cloud is at rest. 
In this frame, the (unnormalized) probability distribution of the states of the single 
particle is 



We can learn various lessons from this. First, there is no Lorentz invariant thermal 
state: a thermal state is in equilibrium in a preferred Lorentz frame of reference, and 
therefore breaks Lorentz invariance. Second, we can say that the form of the statistical 
state is physically determined by the fact that the thermalizing interaction conserves 
p M , and not by the fact that preferred phase space coordinates play an a priori role of 
time and energy. 

The argument above can indeed be sharpened by a more detailed analysis of the 
physics of the system. Note that if the only interaction is elastic scattering among 
the particles, then the gas will diffuse and fail to reach equilibrium. Thus, we need 
something that keeps the gas contained, in order to have a meaningful thermodynamics. 
One possibility to keep the gas contained is to put it in a box. The position of the box 
will then break Lorentz invariance and pick the preferred Lorentz frame in which the 
box does not move. Alternatively, we may think that the particles are gravitationally 
bound. To have Lorentz invariance, we need a field theory for the gravitational field. 
(We can disregard here the difficulties of having point particles, or rigid particles, in 
general relativity, which play no role in this context). In this case, we can approximate 
the dynamics of a single particle as the dynamic of a particle in the mean gravitational 
field of the others. At equilibrium, this gravitational field will be stationary in a a 
preferred Lorentz frame, the one of the center of mass of the cloud. Now, in both cases 
a single particle is not longer free, but rather is subjected to an interaction which is 
not longer Lorentz invariant. In fact, in both cases the interaction preserves energy 
but not momentum (in the second case, energy is kinetic plus potential.) Therefore in 
both cases there is a preferred po which is conserved in the course of the thermalizing 
interaction. It is clearly this energy the one that enters Eq. (|39|) , because this is the 
energy which is totally conserved and freely exchanged in the system, and thus which 
becomes equipartioned. In other words, this is the quantity that is conserved in the 
thermalizing interaction, and that becomes the extensive thermodynamical parameter 
governing the system. 



p(x,p) = e 7MpM 



(37) 




(38) 



p(x,p) = e 



(39) 
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The general lesson should be clear at this point. Even if we do not have an a priori 
recognition of which function on the extended phase space represents time (or energy), 
we can nevertheless run the statistical mechanics formalism on the basis of the quantities 
that are preserved in the thermalizing interaction. 

6. Conclusions and perspectives 

We have argued that quantum statistical mechanical techniques can be applied to 
a macroscopic generally covariant system composed of a large number of generally 
covariant subsystems. This can be done without arbitrarily selecting a variable as the 
time variable, and in spite of the absence of a notion of energy. 

We recall that in the literature there are two main schools of thought in relation to 
the 'problem of time' in generally relativistic theories. One tries to single out the 'correct' 
time variable among the variables of the covariant theory. The choice determines a 
preferred Hamiltonian (energy), and thus an unambiguous concept of temperature. The 
opposite point of view, which we consider more fruitful and we have developed here, takes 
general covariance more seriously, and keeps all variables on the same footing. From 
this point of view, temperature plays no fundamental role in the statistical analysis. It 
may not be defined, or, if it is defined at all, temperature is just one of the intensive 
macroscopic parameters characterizing the equilibrium configuration of the system. 

Our basic idea is that if the macroscopic system can be viewed as being formed 
from weakly interacting systems, then a full solution of the equation of motion of 
the macroscopic system determines a distribution of solutions of the equations of 
motion of the components. The properties of the interaction determine which global 
quantities are conserved and thus the extensive macroscopic parameters. In turn, these 
determine intensive thermodynamical parameters that describe the macrostate. The 
other microscopic degrees of freedom are thermalized away by the interaction. The 
fact that a preferred single notion of temperature does not necessarily arise is not 
surprisingly, given the weak and always contingent role that energy plays in general 
relativistic theories. The statistical mechanics of generally covariant theories does not 
depend on the notion of energy. Rather, ensembles are determined by the properties of 
the thermalizing interaction. 

What is a thermometer in this context? In the usual context a thermometer is a 
physical system which has the property of having a macroscopic variable h (the height 
of the mercury column) directly coupled to the average energy. By looking at h we 
measure directly the average energy and therefore the temperature. If local energy 
is not conserved then a conventional thermometer will keep measuring its own average 
energy, but this will give little information on the system, because individual subsystems 
will not thermalize to the same reading of the thermometer. On the other hand, if other 
intensive parameters 7* are conjugate to other conserved quantities Oi, then in principle 
specific "thermometers" measuring the average value of 0/ may exist. These should 
play the same role as conventional thermometers. 
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Our approach has been very abstract, and applications in realistic general 
relativistic contexts may not be trivial. A first naive idea is to obtain a system 
composed of subsystems by partitioning space into small patches, each with its own 
gravitational degrees of freedom. This procedure, however, might interfere badly with 
general covariance. The example considered in the text suggests to look at the strong 
coupling limit of general relativity, which is precisely given by a collection of finite 
dimensional covariant systems. In this context, we recall that near singularities -such 
as the cosmological one- notions of temperature and entropy are usually very badly 
defined. Alternatively, one could think of somehow Fourier expanding the gravitational 
field, and partitioning it in momentum space, following fluid techniques. Perhaps a 
realistic context in which a covariant statistical mechanics may find application is where 
matter and strong gravitational fields are both present. The presence of matter could 
lead to a natural physical way of partitioning the degrees of freedom. In general, any 
context in which thermal energy can be lost substantially in the gravitational field would, 
in principle, require a covariant thermodynamics. 

On the purely theoretical side, a natural open issue is the quantum statistical 
mechanics of constrained systems with symmetries in the global quantum states. That 
is, the covariant version of Fermi-Dirac, and Bose-Einstein statistics. 

The relation between coherent states in standard quantum mechanics and 
thermodynamics |16| suggests that coherent state quantization of constrained systems 
might bring a better understanding of the thermodynamics of generally covariant 
systems [|T7|], and also shed light on the tantalizing issue of the classical limit states 
in quantum gravity. In particular, consider the spacetime described by a (generic) 
statistical mixture of gravitational states In loop quantum gravity, are 

superpositions of s-knots, or abstract spin networks. In the basis in which p is diagonal, 
we have 

p = £p(V0IVWI (4°) 

and we can compute, for instance, the density matrix yielding average macroscopic 
values of the geometry. In particular, we can use area and volume operators associated 
with a compact region of space and require something like 

A = Tr pA, 

V = Tr pV . (41) 
This may determine a statistical state of the kind 

?=-g-. m 

which might describe the physical state of spacetime better than a somewhat arbitrary 
pure state. 

In closing, let us emphasize again that we have not discussed here the statistical 
mechanics of matter interacting with a fixed gravitational field. Rather, we have 
considered the full quantum statistical mechanics of spacetime itself. Similarly, 
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we have not discussed the statistical mechanics of black holes, which focuses on 
(the perturbations of the gravitational field around) certain preferred black hole 
configuration, and is generally in the asymptotically flat context, in which boundary 
quantities determine time and conserved energy. Finally, we think that the proposal 
of the statistical and algebraic origin of the time flow [|3], |J might be reconsidered at 
the light of the physical ideas discussed here. The magic secret coffer of the relations 
between gravity, quantum and thermodynamics is still far from being fully open. 
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